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We present two results concerning the spin (or pseudo spin) dynamics of 
trapped quantum gases in the hydrodynamic regime described by the Leggett 
equations. First, we apply perturbation theory to extend the "bounded diffu- 
sion" description to trapped systems for small field inhomogeneities. Second, 
we study the formation of long-lived domains with a numerical stability anal- 
ysis of the lowest-lying longitudinal diffusive modes. We also use computer 
simulations ofn/2-pulse experiments to determine the range of experimental 
parameters where the formation of domains can be observed. 



1. INTRODUCTION 

Recent experiments at JILA^ and subsequent analysed have demon- 
strated the identical spin rotation effect (ISRE) in the pseudo spin dynamics 
of trapped ultra-cold atoms with two internal states. The ISRE in spin- 
polarized quantum gaseJ^and analogous mean-field effects in liquid helium 4 
give rise to interesting behavior such as spin waves, instability against trans- 
verse fluctuations, and spontaneous formation of magnetic domains. Fuchs 
et a/P have shown that simulations of the appropriate kinetic equations 
with JILA parameters reproduce the observed pseudo-spin wave oscillations. 
Kuklov and MeyerovictP have noted that Castaing instabilites^ and domain 
wall solutions should be present in trapped g clSGS clS well. Recently, Fuchs 
et a/P have investigated these instabilities by performing simulations of the 
kinetic equations with large initial longitudinal polarization gradients. Al- 
though the Castaing instability is washed out for low collision frequencies, 
they found that the stability should be observeable in trapped gases in, or 
near, the hydrodynamic regime. 

In this paper, we restrict our attention to the hydrodynamic regime 
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(collision time r — > 0) where the equations of motion originally derived by 
Leggett and Ric<P(LR) for spin-polarized systems apply, and where the non- 
linear dynamics are most evident. In the first part, we analyze the behavior 
of trapped systems with small linear and quadratic field inhomogenities, al- 
though the results are readily generalized for arbitrary field inhomogeneities 
(and for finite r, for that matter). In the second part, we carry out a numer- 
ical stability analysis of the longitudinal diffusion modes, and use computer 
simulations to study the formation of domain walls in JILA-type 7r/2-pulse 
experiments. 

For cigar-shaped traps situated along the z-axis the LR equations of 
motion reduce, in the Larmor frame of the center of the trap, to 

/= " ( D fr + nzM ) 4 " d tt^ (i + " Ml x i) < 2 » 

where M = Me, J is the current density, is the linear (k = 1) or quadratic 
(k = 2) strength of the field inhomogeneity, D is the spin diffusion coefficient, 
\iM is the spin-rotation parameter, and = uj 2 t, where uj z is the axial trap 
frequency. We do not take any polarization or density dependence of the 
transport coefficients into account. 

These equations can be put into dimensionless form by writing $7i — > t, 
z^/n/D -»■ z,m= M/Mq, j = J/(M VDn), a = fiM , gi = Gx^DjW, 
and (72 = DG2/Q 2 , where Mq is the equilibrium magnetization at the center 
of the trap. In spherical coordinates, we have 

rh = me + m66 + m sin 0(j>(f) = gkz k m sin 9(f> — j' (3) 



j = — (m + zm) e 



1 + a 2 m 2 



am sin 1 



+ (sin 96 + am6) 



(4) 

where dots (primes) denote temporal (spatial) partial derivatives. After a 
7r/2 pulse the initial conditions are m = exp(— z 2 /2), 8 = tt/2, and 6 = 0. 



2. SMALL FIELDS (ag k < 1) 

For small values of the parameter agk the system is in the bounded 
diffusion regime^ and the external field can be treated as a perturbation of 
the tipped equilibrium magnetization. In this case the field inhomogeneity 
gives rise to (1) transients after the 7r/2-pulse that oscillate at the zero- 
field spin-wave frequencies and (2) a slow decay of the magnetization. The 
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Fig. 1. (a-b) The real and imaginary parts of the eigenvalue spectrum of 
EqEI (c) n = 1 and n = 2 eigenfunctions imp n for a — > oo, (d) Evolution of 
m z at the center of the trap, using only the n=2 mode, for JILA parameters 
n = 1.8 X 10 18 cm~ 3 , u di ff = 0.07Hz (a = 4, g 2 = 0.25). A bias Larmor 
frequency of —2g2 (— 1.5Hz) has been added to simulate a Gaussian field 
inhomogeneity. 



linearized equation of motion is obtained by substituting 9 = tt/2 + 59 (t), 
<p = 5<f)(t) into Eqs.(jHHH), giving 

; , f m V'' \' ■ k i^" ■ k , r > 

mip = Xmip = + igkZ m — > h ig^z m, as a — > oo (5) 

V 1 — iam J a 

where we have introduced the complex phase ip = 59 (t) + i5(p(t). The spec- 
trum of the zero-field modes, shown in Figs, la-b, is obtained by setting 
gk = and solving Eq.((5J) with the boundary conditions imp' — * as z — > oo, 
where m = exp(— z 2 /2). For a = the solutions are diffusive modes of 
alternating parity with A n = —n. For larger values of a some modes of 
opposite parity coalesce into modes localized on the left or right of the trap. 
These modes decay rapidly, however, and are not excited by small fields. 
Only the first two modes with eigenvalues Ai ~ — 1.71/a 2 — 1.32i/a and 
A2 ~ — 7.65/a 2 — A.20i/a are appreciably excited. The analysis is greatly 
simplified for large a in which case the RHS of Eq. (jSJ) is Sturm-Liouville op- 
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erator with kernel m. To analyze the transients after a tt/2 pulse, we write 
tp = Cjity'tfij an d dot Eq.© with ip n , giving 

C'n = iuJ n C n + igkhn (6) 

where bk n = J z k m^ n dz. This has the solution 

c n (t) = ag k d kn (exp iuj n t - 1) (7) 

where dk n = bk n /(&u n ). This solution is still valid for small damping with 
the replacement iuj n — > X n . After the transients the long-lived mode can be 
found from ip = agk^nd-knipn- Alternatively, we set j e = and jg = in 
Eq.(@J). This gives 

9' = am sin 9(f)' = am<fi (8) 

Substituting this into Eqs. (|51Hl gives equations of motion for the long-lived 
mode that are identical to the a = case. These can be solved exactly,^ 
and for small gt- give 5(f) = uJkot + gkZ k , where u>io = and W20 = 52 > with m 
decaying at a rate Ao = g\- 

The m z profile is found by integrating Eq .([%]) and using fm z dz = 0. 
For linear and quadratic fields we get, respectively, 

m z = agxAj^eri ( A=) exp (-^z 2 ) (9) 



m z = ag 2 



1 2^ 1 



cxp(-^ 2 ) (10) 



which are essentially identical to the n = 1 and n = 2 zero-field modes (see 
Fig. lc). We note that for small gt and large a the largest corrections to 
the hydrodynamic results are 0{uj n T oc £lra). Fig. Id shows the evolution 
m z at the center of the trap, which is qualitatively similiar to the evolution 
for larger fields (compare to Fig. 1 in Ref. El where ag 2 = 4.8) . 



3. LARGE FIELDS (ag k > 1) 

For larger values of agk the the modes are strongly coupled and higher- 
orders are evident in the transients. Even so, a linear analysis (which we 
do not provide here) involving the energy of the undamped modes correctly 
predicts n m ax ~ \fotg~k f° r the number of orders observed in simulations (see 
Fig. 2a). Also, for large agk we can treat the trap as a perturbation for 
short times, giving t ~ (a/gk) 1 ^ 2 as the time for the boundary conditions 
to propagate to the center of the trap (i.e. \z\ <C 1). Before this time the 
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Fig. 2. (a) Evolution of m z (z,t) for a = 15 and g\ = 5 showing the higher- 
order transients and the formation of a domain wall, (b) — a stability 
diagram. The small symbols show the results of the numerical stability anal- 
ysis of the n = 1 and n = 2 longitudinal modes. The larger symbols indicate 
whether a domain wall forms (black), or not (white), (gray=marginal) , in 
simulations of 7r/2-pulse experiments. 



magnetization at the center of the trap is effectively unbounded, and for 
large a evolves according to (j> = gizt, 6 = ir/2 (k = 1); or = g2Z 2 t, 
9 = ir/2 + 4g?;t 2 /a (k = 2). 

For large values of ag^ a tt/2 pulse can result in the formation of a 
long-lived domains. These structures, which have been extensively studied in 
liquid helium, 12 ' 13 are formed as the z-component of the spin current j z drives 
the dipole (k=l) or quadrapole (k=2) longitudinal modes. As the transients 
decay, the transverse magnetization {m + = m x + im y ) becomes confined 
near the zeroes of m z . Normally, the field gradient dephases m + rapidly 
and m z decays via ordinary diffusion. However, if a is large enough, the 
longitudinal spin current is unstable, and m + develops into coherent domain 
walls, dividing m z into two (k=l) or three (k=2) domains that decay via 
transverse diffusion across the domain walls. For large ag^ the domain wall 
solution is given b}^ 9" = a(uJk — guz k ) sin#, where u>i = and u>2 = 52^0 ' 
and zq is a zero of the z-magnetization given by erf(zrj/\/2) = 1/2. This 
description is valid so long as the width of the domain wall ~ (agk) 1 ^ 3 "C 1. 

For large gradients the walls are destroyed for g^ > a 3 as found in 
Ref. O for uniform j z . To find the stability threshold for smaller fields, 
we performed a numerical stability analysis of the n = 1 {m z = zm) and 
n = 2 (m z = (1 — z 2 )m) longitudinal modes against transverse perturba- 
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tions. There is a complication that the longitudinal modes decay at rates 
A n = —n. For a given g^, we defined the threshold of stability of the nth 
mode as the value of a where the perturbations grow at the rate +n. The 
results are summarized in Fig. 2b. Determining the minimum a for the 
formation of domain wall after a ir/2 pulse is difficult because the magneti- 
zation decays in a rather complicated manner before the longitudinal modes 
are established. Thus, we carried out simulations of ir/2 pulse experiments 
for various values of a and gt- The presence of a domain wall is signaled 
by its slow decay rate^ A « a -2 / m~ 1 (9') 2 dz k, 2.5<? 1 / 3 /a; 5 / 3 that is much 
slower than the longitudinal diffusion rate. The results, which are summa- 
rized in Fig. 2b, indicate that domains will not be observeable in quadratic 
(or Gaussian) fields for a < 20. However, for linear fields (perhaps in optical 
traps) domains should be observable for a ~ 10. 

Recently, domains have been observed^ in normal trapped 87 Rb, but 
since the lifetime of the domains was comparable to their formation time, 
the experiment corresponds to the (n = 2) marginal category in Fig. 2b. 
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